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Abstract. We study the conditions for operation of the 22 GHz ortho-water maser in a dusty medium near 
late-type stars. The main physical processes, such as exchange of energy between dust and gas in the radiation 
field of a star, radiative cooling by water molecules and pumping of water masers are described self-consistently. 
We show that the presence of dust grains of various types (or of one type with size distribution) strongly affects 
the maser action. The pumping mechanism based on the presence of the dust of different optical properties is able 
to explain water masers in the silicate carbon star V778 Cyg. However, the masers in the winds from asymptotic 
giant branch stars require an additional source of heating, for instance due to the dust drift through the gas. 
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1. Introduction 

22 GHz (6 16 523) ortho- water maser emission has been 
detected from many astronomical objects, such as late- 
type stars, star formation regions, and active galactic nu- 
clei. In recent years, much progress in studying maser 
objects has been made due to increasing angular resolu- 
tion and sensitivity of radio-interferometers, which allow 
one to obtain detailed maps of the distribution of maser 
sources. These observations provide important informa- 
tion about the structure, physical conditions and kinemat- 
ics of the gas in these sources. The observational data show 
that in most of them, maser radiation comes from the re- 
gion co- spatial with sources o f the infr ared dust emission 
( Down es et ahlllisil [Yaniam ura et al ]|2000; Danchi ct al. 

and, therefore, maser amplification may take place 
in a dusty environment. 

In asymptotic giant branch stars (AGBs) and red su- 
pergiants (RSGs), masers are mostly observed in stellar 
winds. Water masers observed in some silicate carbon 
stars, e.g. V778 Cyg (,Engela.l994; Enaels fc Leinert 1994) , 
are not associated with the winds, but most probably arise 
in an oxygen-rich disk around a companion star in a bi- 
nary system ( Szczcrba ct al. 2005). 

The first model of water masers in late-type stars was 
proposed by iDegiichil lll977l) . who considered an expand- 
ing, spherically symmetric circumstellar envelope. The line 
photons were assumed to escape from the wind due to the 
high velocity gradients. The gas temperature was not com- 
puted self-consistently, but was assumed to be a power- 



law function of distance from the star. In his model, the 
strongest water masers are predicted along tangential di- 
rections where the coherent length is largest. He used ap- 
proximate formulae for the collisional cross-sections and 
neglected the influence of dust on the maser strength. 



ICooke fc Elitzuil l)l985(l have presented a more elabo- 
rat e mode l emp loying a new set of collisional cross-sections 
by iGreenI l|l98Ql a nd using the wind paramet ers taken 
from the model of iGoldreich fc Scovilld il97(i\ . In this 
model, the acceleration of the wind is provided by the 
radiation pressure on the dust particles, the temperature 
was assumed to be determined by the stellar radiation, 
while the gas is visc ously heated by c ollisions with faster 
moving dust grains. ICooke fc Elitz^ lll98a assume the 
same temperature profile for outflows with very different 
characteristics. They have accounted for the influence of 
dust emission (using a simplified dust absorption coeffi- 
cient) on the excitation of water molecule, but ignored 
the effect of line photon absorption by the dust. The es- 
cape of photons was again assumed to be determined by 
the velocity gradient. 



iHumphr eys et al"l ( 2001 ) also simulated masers in the 
winds in a ISobole^ ( 196( ) approximation, using more 



recent collisional cross-sections provided by iGreen et alJ 
(,.1993) . They obtained, however, suspiciously high gas 
temperatures because of ignoring molecular cooling. The 
influence of dust absorption on the water molecule popu- 
lation was again neglected. 



2 



Babkovskaia & Poutanen: Water masers in late-type stars 



Recent observations by MERLIN and VLA have shown 
that the winds from AGBs and RSGs are strongly inho- 
mogeneous. Water masers arise at distances of 10-60 AU 
from AGBs and h ave a typical size of 2 — 4 AU across 
l)Bains et al.ll2003^ . while in RSGs all sizes are 5 to 10 
times larger. Moreover, the filling factor of the clouds is 
only ^ 0.01. Because of the relatively small sizes, the 
velocity gradient within the clouds is negligible and the 
Sobolev approximation becomes questionable. 

T he aforementioned rnaser wind models JPesruchil 

ll977tlCooke fc Elitzuill985tlHumDhrevs et al.l20nil) oper- 
ate on the lde Jond lll97.'^ mechanism of coUisional pump- 
ing where the breakdown of thermal equilibrium occurs 
because the gas becomes transparent for different tran- 
sitions at different physical depths. However, the pres- 
ence of dust can strongly influence the maser strength 
llGoldreich fc KwanllT97llKege]l 119751: IStrelnitskiil Il977t 



iBokova et al.lll977j) . lGoldreich fc KwanI I^1974^ suggested 
that the hot dust radiation can excite water molecules 
to the vibrational state, and the heat sink could be 
provided by collisions with cooler (than dust) hydrogen 
molecules. The possib ility of inversio n in such a situa- 
tion was questioned bv lDeguchil l ll98J) . Alter natively, the 
cold dust can produc e the necessary inversi on llStrelnitskii 
T977I iBnlgQvaet a,1.lll977t |neguchilll98n lGha,ndraet a,]' 



IWYTlJlinlgova et aLIIIHYYt |IJeguchilim?^lt K ;ha,nf1ra et, alJ 
1984at ICoHison Watsonlll 99,4 IWallin Watsonlll997l 

m 



Yates et all Il997t iBabkovskaia fc Poutanenl I2004L here- 
after BP04). In an optically thick environment, the exci- 
tation temperature takes a value between the dust and 
the gas temperatures depending on the relative role of the 
dust and collisi ons in the destruction of the line photons. 

has pointed out that water ice can pro- 
vide a very effective heat sink due to a strong peak near 
45 /zm in the absorption coefhcient. For other dust types 
the inversion of maser level populations is smaller by just 
20-30 per cent (BP04). 

In the cold dust-hot gas model, an arbitrarily thick 
layer can participate in the maser action provided the gas 
and dust temperatures differ sufhciently. The temperature 
difference between the gas and the dust can appear in late- 
type stars due to the viscous gas heating by rapidly mov- 
ing dust, presence of the dust particles of different types 
and sizes which assure their different temperatures, or 
shock waves. Thus, the dust present in the outflow plays a 
dual role, providing heating of the gas as well as absorbing 
line photons and increasing the maser strength. Because in 
this situation the photon escape is determined by dust ab- 
sorption and the physical size of the ma ser clump, we can 
use the static approximation ijYates et al. 1997: BP04) for 
maser modeling. This model also can be used for masers in 
a slowly rotating disk around a com panion of the silicate 
carbon star V778 Cyg ijSzczerba et al. 200 5*) . 

The purpose of this paper is to construct a self- 
consistent model of water masers in the environment typi- 
cal for clumpy winds from AGBs and RSGs or maser sites 
in the vicinity of silicate carbon stars. We take into ac- 
count here the most important heating/ cooling processes 
for the gas and dust and study in detail the influence 



of various dust types and sizes on the maser action. In 
Section 2, we formulate the problem and give a detailed 
description of calculation method. Section 3 is devoted 
to the results of calculations and discussion of the multi- 
temperature dust pumping mechanism. In section 4, we 
illustrate our model, applying it to the masers near late- 
type stars. A summary is given in Section 5. 



2. Formulation of the problem and calculation 
method 

We investigate the maser effect in medium containing mix- 
ture of the gas (molecular hydrogen and water vapor) and 
dust in the black body radiation field of a late-type star. 
We simultaneously solve the population balance equation 
for the water molecule and the radiative transfer equation 
(RTE) for water spectral lines using an escape probabil- 
ity method. We also account for the saturation effect in 
masing lines. Our approach is similar to that of BP04. 

In addition, we self-consistently compute the dust and 
gas temperature accounting for main heating and cooling 
processes. As we will see below, the dust temperature is 
determined mainly by the radiation field of a star and is 
independent of the processes in the gas phase, therefore it 
can be computed immediately. 

The gas temperature, however, depends on the energy 
exchange with dust and cooling provided by water va- 
por, which depends on the level populations of the water 
molecule. The populations in turn depend on the gas (and 
dust) temperature. Therefore, we implement the iteration 
scheme where we first guess the gas temperature, compute 
populations and the corresponding cooling rate, then cor- 
rect the gas temperature to satisfy the energy balance, 
until convergence is reached. 

We consider a slab geometry (with half-thickness H) 
which is a reasonably good representation of clumps in the 
winds from AGBs and RSGs flattened in the process of 
formation due to shock waves and/or due to different ex- 
pansion veloci ties in radial and tangential directions (see 
arguments in iBains et al.ll2003l) . It also can be used to 
model masers from disks. The slab thickness 2H is var- 
ied from about 10^^ cm to 10^^ cm which is the upper 
limit corr esponding to the to tal size of the masing region 
in RSGs ijRosen et al.lll97'^ . The stellar temperature is 
taken to be T, = 3000 K. We consider the hydrogen con- 
centration 10"^ < Nii2 < 10^^ cm 
water 10"^ < f^^^ = 9Nii,o/Nii 
is water concentration and the factor of 9 comes from the 
ratio of water to hydrogen masses). 



^ and mass fraction of 
< 10-2 (where iVn^o 



2.1. Radiative transfer 

The solution of the RTE is obtained by the escape prob- 
ability method. The mean intensity of the radiation aver- 
aged over the line profile at the optical depth r (measured 
from the upper boundary of the slab) can be written in 
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the form (see e.g. BP04) 
J(t) = {l-p)Sui+p'B 
p = S + {l-5) 
1 



P 



1 



^K2{0,t) + ^K2{A2to-t) 



1 



(1) 



2L2(/5,r)--L2(/3,2ro-r) 



where 2to = / a^dz is the total slab optical depth, 



Sul = {-, 



1) ^ is the line source function (all intensi- 



ties are measured in units 2hi/^i/c^), Vui is the transition 
frequency, Uu and ni are the fractional level populations, 
gu and gi are the statistical weights for the upper and 
lower levels, respectively, B is the dust source function 
(Planck function), /3 — a^/ai^ is the ratio of the absorp- 
tion coefficients due to the dust and the line, 5{(3) is the 
probability per single interaction that the line photon will 
be absorbed by dust, K2 is the probability that the line 
photon will escape from the slab and L2 is the probabil- 
ity that the photon emitted by the dust will escape from 
the slab. We use approximations of (5, and p'^ from 
BP04. In most of our calculations we take r = tq (i.e. slab 
center). The line absorption coefficient is 



SttcAz/d 



(2) 



Here, An^,; = ni/gi-Uu/gu, /^vb/^'uI ^ {2kT/18mpC^)^^'^ 
is the Doppler width, Aui is the Einstein coefficient, T is 
the gas kinetic temperature, rup is the proton mass and 
At,; = c/ I'ul ■ 

2.2. Masing transitions 

Transitions with inverted level populations (masing lines) 
need a different treatment. We neglect the dust influence 
on masing transitions, i.e. assume P = 6 = p'^ — (e.g. 
|/?| ^ 10""* for the 616 623 transition). The mean in- 
tensity of the radiation in the u/-line in the slab center 
is 



J(r,) 



Sul / [1 - exp(- 
"'0 



(3) 



where fi is the cosine of the angle between the direc- 
tion of propagation and the outward normal, — r^cfix, 
Tin = \/t^Tc — af^H is the optical depth in the maser 
line, Tc is the optical depth in the line center, a™ is the 
maser absorption coefficient, x = {v ^ Vui) / ^i^D is the fre- 
quency within the spectral line in thermal Doppler units, 
— TT^^/-^ exp(— x^) is the Doppler profile. In the case 
of inversion, the optical depth is negative and the mean 
intensity grows to infinity. This happens because the path 
length parallel to the boundary of the slab (/i — 0) be- 
comes infinite. In a real astrophysical system such a sin- 
gularity does not appear because the coherent amplifica- 
tion length is limited by the size of the maser source (or 
velocity gradients). If we approximate the geometry by a 
cylinder with half- height H and radius L, Eq. Q can be 



rewritten in the form 



J{Tx) ~ Sul / [1 - exp(-Tj;/^„iin)] d^ 



Sul / [1 - cxp(-T^//i)] d/i, 



(4) 



where /imin = H/L. The first integral on the right hand 
side is /imiijl — exp(— ra;//iinin)], while the second one is 
approximately [1 - /Zmin - M?nin/kc| exp(-Ta;/Atniin)], when 
\Tx/fJ'mia\ ^ 1- Therefore, for large \Tx/finihi\ the mean 
intensity in the maser line is 



J{Tx) ~ Sul[l - Aiminexp(-rj;//Xmin)]- 



(5) 



At large |T2,//iinin|j Jx grows and the masing transition 
starts affecting level populations (the maser starts to satu- 
rate) . The critical value of at which the maser becomes 
saturated depends on the geometry of the maser source, 
and thus on the value of /imin. The case of /^min = 1 cor- 
responds to an almost spherical maser, while /imin ^ 1 
corresponds to a flat configuration, e.g. a shock wave or 
an accretion disk. 

The line-averaged intensity J = J <f>xJ{Tx)'ix in maser 
lines is given by Eq. (^), where K2 should be replaced by 



One can show that 

A*min exp(|Tc|/^ 

n 



K. 



if 



(6) 



(7) 



while for jr, 
for non 



masing 



0, K™ m ust smoothly approach K^^ 
transitions llHollenbach k McKed Il979 



'Babkovsk aia fc Poutanenll2004|) . We propose here a for- 
mula that connects these two asymptotes (it generalizes 
Eq. [18] in BP04 to ^„ 



< 1 



A'' mi 11 GXp 

n 



1.18|rc 



1 + |rc|[7rln(l + 



11/2 



(8) 



The smoothness of K2, when r passes through zero, allows 
us to account for the effect of maser saturation. 

To characterize the maser strength we use the optical 
depth Till = a™i? determined along the axis of the cylin- 
der. This quantity characterizes only local conditions in 
the slab center. With increasing i?, the maser may stop 
operating in the slab center, reducing such computed t^. 
This, however, just means that maser may operate more 
efficiently away from the central plane closer to the slab 
boundary. 

For a cylindrical maser, the actual observed intensity 
is larger along the radius of cylinder than along its axis. To 
estimate the maximum power of unsaturated cylindrical 
masers, one should use the corresponding optical depth: 



(9) 



If the maser is unsaturated (i.e. the maser radiation does 
not affect the level populations), the absorption coefficient 
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in a masing line does not depend on /imin- However, the 
optical thickness along the line of sight and, therefore, 
the maser luminosity do depend on /imim which is deter- 
mined by the geometry of the source. For simplicity of 
calculations, we assume /imin = 1 everywhere except an 
application in Sec. 13.5.11 

2.3. Population balance 

Using the solution of the RTE in the escape probability 
approximation the statistical balance equations take 
the form (see e.g. BP04) 



^ {piuTii + p'i^BiugiAriui 



1>U 



Ad (puinu + plaBuigiAniu 



Ku 



(10) 



= (Cui^^" - Ciuni) , w = 1, 2, . . . , M - 1, 

with the normalizing condition '^^riu = 1- The Einstein 
^-coefficients corresponding to t he rotational transition s 
of the ground state are taken from lChandra et all l)l984b|) . 
Collisional de-excitation rates with hydrogen, Cui, corre- 
spondin g to the rotational t ransitions of the ground state 
are from lGreen while the excitation rates are 

computed using the detailed balance condition. We take 
into account the first 45 rotational levels of the ground vi- 
brational state and 45 rotational levels of the first-e xcited 
vibra tional state (DIG) of ortho- H2O molecules (iTothI 
Il99lh and all possible radiative and collisional transitions 
between them. We recalculate the Einstein coefficients 
of the rotational transitions in the first-excited vibration 
state and of the vibration-rotational transitions A^j^ 
from the Einstein coefficients of the rotational tran sitions 
in the ground vibrational state A^^ l)Deguchilll977l) 



A., 



4 EG 



A., 



UL 

Mo 



,EG 



(11) 



(12) 



where /x' is the transition dipolc moment between the vi- 
brational states and /io is the intrinsic dipolc moment 
([/i'//io]^ = 0.005), i-'^; is the frequency of a rotational 
transition in the ground vibrational state, is the fre- 
quency of a vibration-rotational transition, is the fre- 
quency of a rotational transition in the first-excited vibra- 
tional state. Since the probabilities of the collisional tran- 
sitions between vibrational states and of the colli- 
sional rotational transitions in the first-excited vibrational 
state CB are not known, we assume, following Deguchi 
l|l977t) . that 



/-.EG 



G 



(13) 



where is the probability of the rotational collisional 
transitions in the ground vibrational state, x = 10~^ is a 
fixed parameter. 



The statistical balance equations (|10|l together with 
Eqs. 1^ for escape probabilities are solved by standard 
Newton-Raphson method. Our escape probability method 
is rather stable and allows us to account for the effect of 
saturation in any masing line, which is very difficult (or 
even impossible) to implement in a n accelerated lambda 
iteration method ijYates et al.lll997j) . 

2.4. Dust properties 

We consider the following types of dust 



amorphous water ice llHud gins et jl1ll993^ . 
crystalline water ice ( Berti(^ ^^L||l9 69l) . 
astronomical silicate (ILaor fc Drain j ll993l ). 
circumste llar silicate (^Davi d fc Pegourielll995() . 
graphite (|Laor fc Draine 1993|) . 



amorphous carbon ( Rouleau fc Martinlll99ll) . 



Water ice is called 'ice' hereafter. 

The dust optical properties can be characterized by the 
absorption coefficient per unit mass ifabs.A = 3QA/4pa, 
where Q\ is the absorption efficiency, p is the density 
of the dust material: 1 gcm~^ (water ice), 2.26 gcm~'^ 
(graphite), 1.85 gcm~^ (amorphous carbon), 3.3 gcm~'^ 
(silicates). We present i^abs,A for different dust types and 
grain sizes as a function of the wavelength in Fig. ^ 
Because of the absence of the experimental data on the op- 
tical constants for the amorphous ice at wavelength A < 1 
/um, we use data for the crystalline ice in this wavelength 
range instead. 

We consider either the fixed size of the dust grains or 
a power-law distribution of the form 



dn{a) oc a''' da, a„ 



< a < a„ 



(14) 



with amin = 0.01 A tm, U maM = 1 /im, and 7 — —3.5 
iMathis c t"anil977l: I.Tural Il996). In our calculations we 
assume the standard dust mass-fraction for interstellar 



medium = 10 ^. Since d/d 



pdn(a)/(2mpA^H2), 



where p is the density of the dust material, one can write 
the distribution function for the dust mass-fraction 



d/d (a) /d 



(7 + 4) 

"max "min 



1''+^ da. 



(15) 



The dust absorption coefficient is given by ad = 
-f^absPd, where pd — /d2?Tip7VH2 is the dust density. If dust 
of two different types is present, the dust absorption coef- 
ficient is given by the sum 



ttd 



2mp7VH. [i^^bs/d + ^abs/d] 



(16) 



where and are the mass fractions of different dust 
types (the total dust mass fraction is /d = /d + /d)- The 
dust source function is then 



B = [alB{Ti) + a'^BiT^)] /ad. 



(17) 
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Fig. 1. Dust absorption coefficient per unit mass, as a function of the wavelength in the case of spherical dust grains 
of radius a = 0.01 uui (solid c urves), a — 0.1 /im (dotte d curves) and a = 1 (dashed curves) for (a) amorphous 
water ice (jHudgins et al.**1993l). (b) crystalline water ice llBertie et al]ll969ll. (c) ast ronomical silicate i Laor fc Draind 
Il993tl. (d) circ mnstcUar sihcate I David fc Pegouri3ll995|) . (e) graphite ijLaor fcPrainc 1993^1 . (f) amorphous carbon 
l|Rouleau fcMartinlll99ltl . 



where T| and are the corresponding dust tempera- 
tures. For the grain size distribution (|15|l . the dust absorp- 
tion coefficient and the source function are transformed to 



ad = 2mpNB_2 J iirabs(a) d/d, 
/ K^hs{a)B[Td{a)]dfd 



B 



J iv:abs(a) d/d 



(18) 
(19) 



We assume that the medium is optically thin to the 
absorption by dust, 



Td = 2mp/diVH,i^abs(r,)i/ < 1, 



(20) 



where Td is the dust optical depth, -fCabs is the Planck- 
averaged absorption coefficient. 

2.5. Dust cooling and heating 

The main mechanism for dust heating is provided by the 
stellar radiation 



In addition, there is energy exchange between gas and dust 
(Goldrcich & Kwan 1974) 

<Zgas = NdAdNn,Vtk{T-Td) = ^^Nl,Vtk{T^Td), (23) 

where A^d is the concentration of the dust grains, = ira^ 
is their geometrical cross sections, and Vt is the hydrogen 
thermal velocity. The energy exchange is small compared 
to radiative cooling and heating if A^H2 ^ 10^ cm^'^ and 
can be neglected. Absorption of the near-infrared photons 
from H2O and CO present in the dusty medium gives an 
even smaller contribution to the energy exchange. The 
thermal balance for the dust is then transformed to a stan- 
dard form 



'Jrad- 



(24) 



The dust temperatures depend only on the stellar temper- 
ature, the dilution factor (approximately as VF^/^) and the 
dust optical properties (see Fig. 



(21) 2.6. Gas heating and cooling 



where = i (g) is the dilution factor, R is the stellar The main heating and cooling mechanisms determin- 
radius, and D is the distance to the star. The dust cools ing the gas kinetic temperature are th e dust-gas col- 
by its own radiation with the corresponding rate lisions, photoelectric effect on grains llde Jond |l977t 

iTielens fc HoUenbach 1985; Groenewe genI Il994l) . cos- 
9rad = 8mp7VH2/d-fi^abs(Td)crr^. (22) mic rays (,Goldsmith fc Langer. 119781) . viscous stresses 
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Fig. 2. Dust temperature for dust grains of radius a — 
0.01 /^m for amorphous ice (dotted curve), crystalline ice 
(solid curve) , astronomical silicate (dashed curve) , circum- 
stellar silicate (dot-dashed curve), graphite (triple-dot- 
dashed curve), amorphous carbon (long-dashed curve), as 
a function of the dilution factor for T, ~ 3000 K. 



ijShakura fc Suiivaevlll973l) and radiation o f H2O and H2 
molecules fe.g. iGoldreich fc Scovill^llQTGfl . The heating 
(cooling) rate of the gas per unit volume (ergcm~'^s~-^) 
due to collisions with hotter (cooler) dust is given by 
Eq. (123: 

^dust = "'J'gas- (25) 



The rate of gas heating due to cosmic rays and viscous 
stresses is several orders of magnitude smaller than that 
due to the dust-gas collisions, while heating of the gas due 
to the photo-effect is important when the stellar radiation 
is strong in the UV range, and can be neglected for late- 
type stars. Moreover, in spite of the fact that hydrogen is 
the most abundant molecule in the oxygen-rich material, 
the rate of gas cooling due to H2 (Hartauist et al. 198(|) 
is much smaller than that due to the dust-gas collisions. 

Next to H2 the most abundant molecule in 
oxygen-rich matter is prob ably the H 2O m olecul e 
llGoldreich fc Scovilld Il976t K^roeneweeenl Il994|) . 
iNeufeld et alJ l)l995l) obtain cooling functions for 
five molecules and two atomic species that are potentially 
important coolants (H2, H2O, CO, O2, HCl, C, and O), 
and show that, for a temperature of about 200 K and 
hydrogen concentration 10^ < A^Hs ^ 10^ cm~'^, the gas 
cools mostly by H2O molecules. The heat loss rate of the 
gas due to radiative cooling by H2O is 



9H2O 



(26) 



where the sum is over all collisional transitions in the wa- 
ter molecule, Ciu and Cui are the rates of the collisional ex- 




Fig. 3. Water cooling efhciency (erg cm^ s~^) as a func- 
tion of Nh^Nh^qH in the case of the water- hydrogen mix- 
ture, including the first vibrational level in water molecule: 
with silicate dust (solid curve) and without dust (dashed 
curve); taking into account only ground vibrational state: 
with silicate dust (dotted curve) and without dust (dot- 
dashed curve). The dust and gas temperatures are 200 
K, the grain size is 0.01 /im, dust-to-gas mass ratio is 
/d — 10^^, the slab half-thickness is H = 10^"* cm. 



citation and de-excitation. In principle, the value of gHjO 
can be negative, which means that the radiative cooling 
in water lines is less than the heating due to absorption 
of dust radiation. 

In order to test our code, we compute the cooling 
rate due to water lines. Fig. 13 shows the dependence 
of gH2o/(-^H2-^H2o) on NH2NH20H. First, we take into 
account only the ground vibrational state in the wa- 
ter molecule and then include the first-excited vibra- 
tional level. Our results a re similar to those obtained by 
iNeufeld fc Melnickl l|l987t ). who accounted for the first 179 
rotational levels of the ground vibrational state. The small 
difference can be explained by the different number of 
rotational levels included. In the presence of dust with 
Td = T the cooling is reduced because dust grains trap 
part of the radiation from water molecules in the medium. 
If one includes the first-exited vibrational state in the wa- 
ter molecule, the cooling rate becomes larger because of 
additional vibration-rotational collisional transitions and 
rotational transitions in the first-exited vibrational state, 
but with dust there is no difference. 

A rather general form of the thermal balance equation 
for the gas is 



<Zdust — <ZH2 0- 



(27) 
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The thermal balance for the mixture of the gas and dif- 
ferent dust types (noted by index i) is 



(28) 



The gas temperature can be obtained by simultaneously 
solving this equation with the population balance equation 
for a water molecule (see section IT^ . Such a system of 
equations is non-linear and is solved by iterations. If the 
radiative cooling by H2O is small, the gas temperature is 
expressed by a simple relation 



T = 



J2^ Nd,iAd,t Ei fda/Piai 



(29) 



In the case of the grain size distribution, the energy bal- 
ance for the gas is 



3mp 



d a 



T]a M/d = gHaO- 



(30) 



When various dust types are present, some can cool the 
gas while others can heat it. 

We illustrate the effect of water cooling with a simple 
example of two dust types mixed with hydrogen and wa- 
ter. We consider equal mass fractions of circumstellar sili- 
cate and amorphous carbon with total /d = 10~^. Silicate 
is cooler, while carbon is hotter. It is clear from Eq. (|30|l 
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Fig. 4. Dependence of the ratio '/Hao/'Zdust '^^ dilution 
factor W for a circumstellar silicate and amorphous car- 
bon dust mixture at the fixed values of slab half-thickness 
H = 10^'' cm, the gas concentration A^H2 = 10* cm~^, 
water mass fraction /H2O = 10^^, dust mass fraction 
/d = 10~^, temperature of the star = 3000 K and 
different grain size a — 0.01 /^m (solid curve), 0.1 /^m 
(dotted curve), 1 //m (dashed curve), assuming the same 
mass fraction of the two dust types. 



100 



a=0.01/xm 
I I 




L a=0.1/^m _E 

j I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I r 




a- l/i,m 



10^ 10^ 10'' 10' 

II cm 



10 



10 



AAfj [cm' 



Fig. 5. Gas temperature as a function of hydrogen con- 
centration in a slab of half-thickness H ~ 10^^ cm in 
the radiation field with W = b x 10"^ and = 3000 
K and diff'erent grain sizes (a) a = 0.01 /im, (b) 0.1 /xm, 
(c) 1 fim. Water fraction of /hsO = 10^^ (solid curves) 
and /H2O — 10^'^ (dot-dashed curves) is assumed. Dust of 
two types, silicate and amorphous carbon, is present with 
equal mass fraction and total /d — 10^^. Dust tempera- 
tures of the silicate and amorphous carbon are shown by 
dotted and dashed lines, respectively. 



that small grains are more effective in heating (and cool- 
ing) the gas. Fig. 01 compares the water cooling rate to 
that due to collisions with silicates for different grain sizes 
a = 0.01, 0.1, and 1 /zm as a function of the dilution fac- 
tor. We see that the radiative cooling rate due to H2O is 
about two orders of magnitude smaller than the cooling 
due to collisions with 0.01 /im grains, while the efficiency 
of these two mechanisms is comparable for 1 fim grains. 

The effect of radiative cooling by H2O on the gas tem- 
perature does not depend only on the size of the dust 
grains, but also on the hydrogen concentration. This is 
because the heating and cooling rates have different de- 
pendencies onNn^: g'dust(-^H2) -^Ha/d (seeEqs. ^51 150) 1 

and gH2o(A^H2) « (iVH2iVH2o)'/' « N^J^H (see Fig.E)). 
In Fig. O we illustrate the dependence of T on N^^ for 
the dilution factor = 5 x 10~^ (where water cooling is 
most effective, see Fig.^J. At high Nb_2 the gas tempera- 
ture takes an intermediate value given by Eq. (|29|l . In the 
case of 0.01 /im grains, H2O radiative cooling is negligi- 
ble for Nb_2 ^ 10^ cm~^, while in the case of 0.1 and 1 
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fim grains, water cooling is important even for iVua ^ 10^ 
cni~^. Moreover, for large grains and A^H2 ^ 10^ cm~'^ 
water cooling is so effective that the gas becomes cooler 
than the cold dust (silicate). 

In the winds from late-type stars, the radiation pres- 
sure from a star accelerates the dust which can drift 
through the gas, providi ng an additional source of gas 
heating (see, for example. iGoldreich fc Scovillelll976() 

1 

2^ 



Adrift ^ ^Nii^2mpYNd,iAd,iVli, 



(31) 



where Vd is the drift velocity. This term should be added 
to the balance equations lP7|l and 

3. Results 

3.1. de Jong model 

Let us first discuss how the 6i6 — 623 water maser works 
if no dust is present in the medium. We will thus have 
a benchmark to which we can compare our dusty mod- 
els. De Jong (1973) considered a gas cloud consisting of a 
mixture of water vapor and hydrogen. He showed that 
upon approach to the surface of the cloud the optical 
depth in the lines decreases and the thermal equilibrium 
breaks down. This occurs at different depths in differ- 
ent lines and over some range of depths some levels be- 
come relatively overpopulated, leading to population in- 
version in some transitions. The coUisional rates of the wa- 
ter molecule with hydrogen were not known at that time, 
and he used approximate expressions. We rep eat these cal- 
culatio ns using modern coUisional rates from lGreen et al] 
l|l993(l . To find the main transitions participating in the 
maser pumping mechanism we compared the population 
fluxes Fui = UuWui - niWiu (where Wui and Wiu are to- 
tal transition probabilities) for all possible pure rotational 
and ro- vibrational transitions of ortho-H20 molecule in 
the ground and first-excited vibrational state. 

We find that for N^^ > 10^ cm'^ and /h^o > 10"^ 

the mam cycle ot the maser pumping is 023 — > 032 — > 
641 ^iiilf 616 (solid lines in Fig.l^. The upward transitions 
523 532 and 832 641 are dominated by collisions and 
therefore, water molecules are coUisionally excited from 
523 to 541 level through 832 level. Radiative de-excitation 
from 541 to 6 16 produces the heat sink in the maser pump- 
ing cycle, when radiative rate of this transition is larger 
than the coUisional one and the optical depth in the cor- 
responding line is small. In the remaining range of param- 
eters iVHa and /H2OJ the main cycle of the maser pumping 

IS 023 — > 4i4 — > 7o7 — > Dig. Water molecules are ra- 
diatively de-excited from the 823 to 4i4 level, coUisionally 
excited to the 7o7 level, and again radiatively de-excited to 
616 closing the pumping cycle. The pumping cycles found 
by us appear to be different from those considered by de 
Jong, probably because of more accurate coUisional rates 
in our calculations. 

In Fig. [3 (case 1, dashed curves) we present the 
maser absorption coefficient as a function of the slab 
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Fig. 7. The optical depth in the 616 — 623 maser line as 
a function of the half-thickness of the slab for Nb_2 = 10* 
cm~'^, /H2O = 10"'^. The following cases are considered: 
(1) de Jong mechanism in the water- hydrogen mixture 
(without dust) of T = 230 K (dashed curves); (2) same as 
previous, but with circumstellar silicate grains of a = 0.01 
/im and = T (dotted curves); (3) gas mixture with 
amorphous carbon grains of a = 1 /im and = 279 K. 
The gas temperature is calculated self-consistently, using 
thermal balance equation H28|l (dot-dashed curves); (4) 
circumstellar silicate grains of a — 0.01 /im and = 143 
K, gas with T = 230 K, and 0.01/im amorphous carbon 
grains with = 279 K (solid curves). Dust temperatures 
correspond to the dilution factor W = 5x 10"^ and stellar 
temperature = 3000 K. We assume = /| = 
/d = 10~^, and /i,„in — 1- On the right hand side of the 
vertical dashed lines the dust is optically thick. 



half-thickness. The maser disappears at H > 10^^ cm, if 
A^H2 = 10* cm~^, when the probability of photon escape 
in 541 — > 616 transition (heat sink) becomes comparable 
to the rate of collisions, i.e. Cs^^^gie — P^54i-»6i6 (see 
Eqs. |lllOj and BP04). For the same water concentration 
but higher hydrogen density of Nb_2 — 10^ cm~^, this 
condition is satisfied at smaller line optical depth r corre- 
sponding to H ^ 10^^ cm, where the maser disappears. 

3.2. Influence of dust 

Let us now investigate the infiuence of the dust on the 
inversion of the maser level populations. If dust (we con- 
sider 0.01 /im circumstellar silicate grains) is of the same 
temperature as the gas, the maser efficiency drops at large 
depths (compare case 2, dotted curves in Fig. |7|to case 1, 
dashed curves). This happens because dust traps radia- 
tion from water molecules in the medium, while the de 
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Fig. 6. System of the rotational energy levels of ortlio-H20 molecule in the ground and first-excited (010) vibrational 
state. Arrows show the main cycles of the maser pumping in the following cases: Deguchi model (dot-dashed lines), de 
Jong model (solid lines), two-dust model (dashed lines). The cycle extinguishing the maser effect at large H for the 
combinations of the ice with any other dust is shown by dotted lines. 



Jong mechanism is based on the escape of photons from 
the surface. 

We have not specified any mechanism of gas heating. 
In the models considered above we have assumed certain 
gas and dust temperatures without considering energy bal- 
ance. If there is no relative motion between the dust and 
the gas, the only efficient mechanism of gas heating is 
the energy exchange with the dust. If only one type of 
dust is present, the gas temperature can be at most that 
of the dust, and if water cooling is important then it is 
much smaller. The hotter dust not only heats the gas, but 
strongly affects the maser operation because of the influ- 
en ces of the dust radiation. 

iGoldreich KwanI l(T97i proposed that the inversion 
of the maser level populations 6i6-523 can arise due to the 
excitation of the first-excited vibration state (010) of the 
ortho-H20 molecule by 6.3 fim dust radiation, while the 
heat sink is realized through collisions w ith coo ler hydro- 
gen. This model has been criticized by IPeguchl , 



who has pointed out that at high concentrations needed 
for the coUisional sink, the levels at the ground vibrational 
state will be thermalized. 

To check this scenario, we considered amorphous car- 
bon grains of size a = 1 fim in the radiation field with 
stellar temperature = 3000 K and dilution factor 
W = 5 X IQ-^. The dust temperature Td = 279 K is 
determined from the thermal balance equation H24|l . The 
gas temperature, which is computed self-consistently from 
Eq. 128() accounting for the water cooling and the heating 
by collisions with the dust, varies from 43 K at iJ = 10^^ 
cm to 276 K at iJ = lO^^ cm (for N^^ = 10** cm^^). The 
inversion of maser levels (see case 3, dot-dashed curves in 
Fig. O is smaller than in previously considered cases and 
is still produced by the de Jong mechanism. At small H 
the radiative cooling by water is so intense that the gas 
temperature becomes too low for maser pumping. In spite 
of inability of the hot dust to efficiently pump the maser. 
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Fig. 8. The optical depth in the 6i6 — > 623 transition 
as a function of the ratio //j for different cold - hot dust 
combinations: amorphous ice and graphite (1), crystalline 
ice and graphite (2), crystalline ice and amorphous carbon 
(3), crystalline ice and circumstellar silicate (4), circum- 
stellar silicate and graphite (5), astronomical silicate and 
graphite (6), circumstellar silicate and amorphous carbon 
(7), astronomical silicate and amorphous carbon (8). The 
parameters are H = lO" cm, = 5 x 10"^, = 3000 
K and the grain size a = 0.01 /xm. 



log /h,0 

Fig. 9. Contour plots of levels of constant the 616 623 
maser optical depth —t^. Solid curves correspond to the 
case of the circumstellar silicate dust with = 143 K 
and amorphous carbon with = 279 K; dotted curves 
correspond to the case of equal gas and dust temperatures 
T = Td = 230 K. The size of the dust grains is a = 0.01 
lim and the slab half thickness is H = lO^'*^ cm, the mass 
fraction of the silicate and carbon dust is the same = 
= 5 X 10""^. On the upper side of the horizontal dashed 
lines the dust is optically thick. 



the presence of the hot dust settles the problem of the gas 
heating. 

Of course, if there are additional mechanisms of gas 
heating (e.g. due to dynamic friction with rapidly mov- 
ing dust, see Eq. the gas temperature can exceed 
the dust temperature by a large factor. In this situa- 
tion the maser is very effective (as was discussed for ex- 
ample by 'Chan dra et alJll984at IColHson &: WatsonlllQQSt 
lYates et al.lll997t BPOA). 

3.3. Hot and cold dust pumping mechanism 

The lat e-type stars often sh ow emission of different types 
of dust ((Danchi et alJll994l) . Thus, it is natural to assume 
that two dust types are present simultaneously. Since dust 
temperatures differ (see Fig. collisions with one dust 
can heat the gas, while with another can cool it. The gas 
temperature then takes an intermediate value (see Fig. [S]) . 

What pairs of dust produce the largest inversion? From 
Fig. [3 we see that a significant temperature difference can 
be achieved for dust types where one is graphite or amor- 
phous carbon and another one is any of the remaining (sil- 
icates or ice). For completeness, we also consider a mixture 
of crystalline ice and circumstellar silicate, in spite of the 
fact that the temperature difference is not so large. Fig.|Sl 



shows the dependence of the optical depth in the maser 
line on the ratio of the hot-to-cold mass fractions //d 
for different dust combinations. We see that the strongest 
maser (i.e. largest negative optical depth Tm) occurs at a 
ratio /d//d — 1 foi' combinations involving graphite (or 
amorphous carbon) and silicates. We further assume the 
same mass fraction of the cold and hot dust, /d = /d- 

Let us now consider a mixture of 0.01 /im circumstel- 
lar silicate and amorphous carbon grains. Their temper- 
atures, in the radiation field with W = 5 x 10^^ and 

= 3000 K, are T^^ = 143 K and T^" = 279 K, re- 
spectively. The gas temperature of T = 230 K is fully 
determined by the collisions with dust which arc more ef- 
ficient than water radiative cooling for small grains (see 
section The resulting maser strength as a function of 
the slab thickness is shown in Fig. [T] (case 4, solid curves). 
We see that aX H < 10^** cm, the inversion is identical to 
that in the dust-free case, implying operation of the de 
Jong mechanism. At larger H, the maser absorption coef- 
ficient is almost constant resulting in a linear increase of 
the maser optical depth. At even higher H >3 x 10^^ cm 
the maser saturates (we used fimin = 1 in these calcula- 
tions, i.e. almost spherical maser, see Eq.lSjl. Thus, we see 
that the presence of two dust types is rather efficient in 
heating the gas as well as keeping a large inversion deep in 
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Fig. 10. The optical depth in the 6i6 — 623 maser hne 
as a function of the slab half-thickness for the same cold- 
hot dust combination as in Fig. |S1 and the following set of 
parameters: a = 0.01 /xm, = 3000 K, VK = 5 x 10-^ 
/d = /d = 5 X 10-3. 



the cloud. For N-n^H > 10^"^ cm^^ the medium becomes 
optically thick to absorption by dust, and a large differ- 
ence between the gas and dust temperatures is unlikely 
to be produced in real astrophysical environments (unless 
there is additional heating, of course). 

In Fig. Owe present contour plots of levels of the con- 
stant 616 — > 523 maser optical depth T,n at the plane hy- 
drogen concentration ~ water-to-gas mass ratio. The gas 
temperature is calculated self-consistently from the ther- 
mal balance equation H28|) . For large N-^^ or small water 
content, the gas temperature is 230 K, while for smaller 
TVhj and large /H2O, the gas cooling by H2O is impor- 
tant and T decreases down to 150 K (lower right corner in 
Fig. 1^1). We assume that the dust temperatures are con- 
stant in the whole range of A^Ha • However, for TVhs ?J 10^ 
cm^"^ (and our assumed H — 10^^ cm) the medium be- 
comes optically thick to the dust absorption and therefore 
the gas and dust temperatures equalize. To investigate 
the maser effect in the case of equal gas and dust tem- 
peratures, we also present in Fig. |^ (dotted curves) the 
maser optical depth for T = = 230 K. In that case the 
inversion of maser levels exists up to N-^^ < 10^" cm"'^ 
and the maser is weak in the range 10^ < N-^^ ^ 10^" 
cm~3. Thus, the 616 — > 523 maser is mostly effective in 
the medium with < 1 (see Eq. |20j). 

The 6 16 — 523 maser optical depth as a function of H is 
shown in Fig. ^| for the same cold- hot dust combinations 
as in Fig. |S| We see that the maser is less effective for the 
combinations of ice with other dust types and most effec- 
tive for the combinations of silicates with amorphous car- 
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Fig. 11. The 616 — 523 maser optical depth as a func- 
tion of the dilution factor W for different cold-hot dust 
combinations: astronomical silicate and amorphous car- 
bon (solid curves) , astronomical silicate and graphite (dot- 
ted curves), circumstellar silicate and amorphous car- 
bon (dashed curves), circumstellar silicate and graphite 
(dot-dashed curves), crystalline ice and amorphous car- 
bon (triple-dot-dashed curves), crystalline ice and circum- 
stellar silicate (long-dashed curves). The parameters are 



/d = /d = 5 X 10" 



0.01 fim, H = lO^'' cm, and 



T* = 3000 K. Maser efficiency for ice is not computed for 
W > 10"'*, because of its evaporation at corresponding 
temperatures. 



bon. The maser works due to the de Jong mechanism when 
H < lO*"* cm, while for the silicates- carbon (graphite) 
mixtures the inversion exists also at larger H, where the 
maser action is controlled by the dust. 

The analysis of the population fluxes F^i reveals that, 
unlike the cold dust - hot gas mechanism operating on 

the pumpmg c ycle 523 — > 4i4 — > 5o5 — > Oie 
l|Deguchil Il98lt BP04; see dot-dashed lines in Fig. EJ, 
the main pumping cycle of the two-dust mechanism is 

-"34 ''^^^ 743 ^-^^ 616 (dashed lines in Fig. EJ. 



J23 



5.5;jm g 
> Do 



In other words, the presence of the hot dust extinguishes 
maser pumping based on the absorption of 45 fim pho- 
tons, because the absorption coefficients of silicates and 
amorphous carbon near 45 /zm are similar (see Fig.0. 

Since the absorption coefRcient of the amorphous car- 
bon at 5.5 /Ltm is much larger than that of silicates, the 
523 ^^^^ transition is dominated by the hot dust. On 




Fig. 12. (a) Dependence of the dust temperature on the grain size for graphite (sohd curve), amorphous carbon (dotted 
curve), circumsteUar siUcate (dashed curve), astronomical sihcate (dot-dashed curve), crystaUine ice (triple-dot-dashed 
curve), amorphous ice (long dashed curve), (b) Gas temperature as a function of the hydrogen concentration for grains 
of the same dust types as in (a) with the size distribution described by Eq. Other parameters are /naO = 10"'^, 
/d = 10-2, H = 10" cm, = 5 X 10"^ T* = 3000 K. 



the other hand, there is a peak in the absorption coef- 
ficient of the silicate dust near 20 fim (see Fig. and, 

therefore, the 743 transition is dominated by the 

cold dust. Thus, level is populated by absorption of 
the 5.5 /im radiation from the hot dust, while the heat 
sink is realized by absorbing the 21 /im photons by the 
cold dust. Moreover, since the Einstein coefficient for the 
634 ~* ^43 transition is large A — 2.4 s~^, water molecules 
are radiatively de-excited from Gf^ to level 743, closing the 
cycle of 616 523 maser pumping. 

When cold ice is combined with any other hot dust the 
maser effect is extinguished at large H because of the cy- 

, ^ 30uni _ 75um ^ 49um _ / i , , it . x^. im 

cle 616 — > 125 — > 634 — > 523 (dotted hues m l<ig.|6|). 
This cycle arises due to the gap in the ice absorption co- 
efficient near 30 fj,m (see Fig. The molecules are ex- 
cited from 616 to 725 by absorption of 30 /im photons 
from the hot dust, which dominates this transition. De- 
excitation from 725 to 523 is realized through 725 ^-^^^ 634 

and 634 transitions dominated by the cold ice 

(since its absorption coefhcient is larger than that of sili- 
cates or amorphous carbon at these wavelengths) and in- 
version disappears. 

Fig. shows the dependence of the maser power on 
the dilution factor for the dust combinations (3)-(8) con- 
sidered in Figs. IHl El One can see that the amorphous 
carbon appears to be more effective than graphite in pro- 
ducing the inversion of the maser levels, while the astro- 
nomical silicate is more effective than the circumsteUar 
silicate. It happens because the absorption coefficient at 



21 /im (743 — > 616 transition responsible for the heat 
sink in the maser pumping cycle) of the amorphous car- 
bon is smaller than that of graphite, while the absorption 
coefficient of astronomical silicate is larger than that of 
circumsteUar silicate (see Fig.^. The larger the efficiency 
of absorption of these photons by the cold dust and the 
weaker the emission of these photons by the hot dust, the 
larger the inversion of the maser level populations. 



3.4. Grain size distribution 

In the previous section we showed that the presence of two 
types of dust solves the problem of gas heating as well as 
provides the necessary conditions for strong maser pump- 
ing. A similar effect can be achieved if there is a dust grain 
size distribution. In Fig. 112b we present the dependence of 
the dust temperature on grain size. We see that the sili- 
cate, graphite and amorphous carbon temperatures vary a 
lot, while the ice temperature is almost constant, because 
optical properties of the ice weakly depend on the grain 
size (see Fig. ^ . The radiative cooling by water molecules 
strongly affects the gas temperature at small N^i^ where 
collisions with dust are less efficient (see Fig. 112b'). At 
A^Ha ^ 10* cm""^ the gas temperature approaches the dust 
temperature of the smallest grains which dominate the gas 
heating (see Eg. I30|l . At N-^^ ~ 10^ cm^'^ (for silicates at 
10* cm"'^), the gas temperature becomes smaller than the 
minimum dust temperature. 
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Fig. 13. The optical depth in the 6i6 — 623 maser hne as a 
function of the half-thickness of the slab for the same dust 
types as in Fig. 1121 and the same parameters, Nn^ = 10^ 
cm~^. 
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Fig. 14. Contour plots of levels of constant 616 — > 623 
maser optical depth, — t^, for the amorphous carbon dust 
distribution (|14|l . Same parameters as in Fig. 1121 



temperature thus becomes too low for maser pumping and 
the maser disappears. For ice, silicate and amorphous car- 
bon the inversion exists only for 10^^ H < 10^^ cm, 
which means that here the maser is pumped due to the 
de Jong mechanism. Graphite and carbon dust have the 
highest temperatures and largest temperature variations. 
This causes the maser to operate even in the deep layers. 
We see also that the maser is stronger for ice dust than for 
the silicate, in spite of the fact that the ice temperature 
is smaller than of the silicates (see Fig. \12h). The reason 
is that the water cooling rate is larger in the presence of 
silicates than in the presence of ice. It happens because 
the ice traps radiation better than silicates, since H2O 
lines are mostly concentrated in the IR region, where the 
absorption coefficient of ice is larger than that of silicates. 

The dependence of maser optical depth on the hy- 
drogen concentration and water-to-gas mass ratio is shown 
in Fig. ll4l for the amorphous carbon dust distribution. The 
maser strength is large near N^^ — 10^ -f- 10^ cm^'^ and 
high water concentration. These results arc similar to that 
for the case of the two-dust model (see Fig. because 
the cycle of the maser pumping in these two models is the 
same. 



3.5. Applications 

3.5.1. Silicate carbon star V778 Cyg 

Water masers are observed in some silicate carbon stars, 
e.g. V778 Cyg llEngelalll994l:lEngels fc LeinertJll994l) . This 
object also shows silicate dust em ission features in the 
IR spectra l|Yamamura et al.ll200r)|) . Recent detailed map- 
ping of the masers with MERLIN reve aled an elongate d 
S-shape structure of about 32 AU ( Szczerba et alJl200!T|) . 
Such a structure can be interpreted as a warped disk 
around a low-mass companion of the carbon star in a bi- 
nary system observed almost edge-on. 

The lower limit on the maser brightness temperature 
is Tb = 6 X 10* K. This requires the maser optical depth 
of at least t — — ln[rR/T nx] ^ —15.6, where the excitation 
temperature is ~ 100 K. ISzczerba et al.l l)2005|) estimated 
the distance between the masing disk and the carbon star 
to be about 80 AU. With the stellar radius of 2 AU, the 
dilution factor is about W — 10~^. 

We assume that the disk contains a mixture of O.Ol/im 
astronomical silicate and amorphous carbon grains with 
the total mass fraction /d = 0.01. In the stellar radiation 
field of a carbon star (with temperature of about 3000 K), 
the temperature of the dust is Tas = 235 K and Tac = 500 



In Fig. ^1 we present the dependence of the maser 
optical depth as a function of the slab half-thickness. 
The maser absorption coefficient is maximal for H = 



10 cm and water frac- 



K, respectively. Taking 

tion /H2O = 6 X 10"'^ IjJeong et a l.''2003V the resulting gas 
temperature is T = 360 K. From Fig. 1101 we get t,„ « 8 
for H = 10^^ cm which gives us the absorption coefficient 
a™ = T^/H = — 8 X lO^^'' cm^^. For the projected disk 
10^^ cm and disappears at H < 10^^ cm. It can radius of i? = 16 AU, the coherent length in the disk is 



10' 

be explained as follows. The smaller H , the smaller is the 
trapping of photons in the slab and the larger is the ra- 
diative cooling of the gas by water. At small H the gas 



about S = 1.2R = 3 x 10^"' cm. Using Eq. we get 



= S 



-13.5. However, for H — 5 x 10^^ cm 



(e.g. away from the disk central plane), the inversion is 
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larger by a factor 1.5 (see Fig. Ill|l . This then gives the 
maximum maser optical depth of —20 which is larger than 
the observed lower limit. Even for a much lower water con- 
tent /H2O = 10"'^ we still get r = — 14. Thus, coUisional 
energy exchange with the dust can provide the gas heating 
which is necessary for the masers to operate. 



3.5.2. Masers from AGB winds 

Water masers from AGB stars are observed in their ex- 
panding envelopes. Four low-mass late-type stars (IKTau, 
U Ori, R TVir and U Her) have been mapped by MERLIN 
recently l|Bains et al.1l200^ . The data show that maser ra- 
diation comes from individual maser clouds with the ap- 
parent size of 2 — 4 AU (Richards et al, 1999) and the 
filling factor of only ^ 0.01. Masers are observed at typi- 
cal distance from the star of 10 — 70 AU. The total maser 
photon production rate (luminosity) is about (1 — 8) x 10^^ 
s~^. With the total number of clouds varying between 14 
and 286 depending on a source, one can estimate a sin- 
gle cloud luminosity of less than 10^"'^ s^^. For a spherical 
cloud with the radius of ~ 2.5 AU, the maser emissivity 
$ varies between ~ 0.2 and ^ 0.6 cm^'^ s"^. 

Let us consider a gas - dust cloud of radius H — Ax 10^^ 
cm in a radiation field of the star with = 3000 K and 
W = 10~4 (i.e. at a distance of about 50 AU for a 1 
AU stellar radius), and other parameters are the same 
as in section r^. 5. II The gas temperature depends on the 
heating which can be provided by collisions with dust. We 
consider heating (a) due to the gas thermal motion (see 
Eq-ESI) and (b) due to the drift of grains through the gas 
(see Eq.PH). If only heating due to thermal motions is 
considered, the gas temperature is about 360 K and it lies 
between temperatures of the cold and hot dust. The maser 
optical depth is Tm ~ —5 (see dotted curve in Fig. I15|) . 
From our model we also c an compute t he average maser 
photon emissivity (see e.g. lElitzuilll99l[l 
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Fig. 15. Dependence of the maser optical depth, maser 
photon emissivity $ [photon cm~^s~^] and temperatures 
of the gas and dust on a height z within a slab (z = in the 
slab center). Solid curves show the case of the drift heating 
of the gas given by Eq. and dotted curves represent 
the case of heating due to thermal motions (see Eq. [23 )■ 
Dashed and dot-dashed curves represent the temperatures 
of the astronomical silicate and amorphous carbon, re- 
spectively, calculated for W = lO"'^ and = 3000 K. 
Other parameters are 7Vh2 — 10* cm^'^, /H2O = 6 x 10"'^, 
/d = 10-2, H = Ax 1013 cm, a = 0.01 /xm, drift velocity 
Vd = 2 km s-i. 



$ = guA^i \An^i\ Nh,o\Sui\ K^, (32) 

which is about 2 x 10~"^ cm"'^ s-^, i.e. two orders of mag- 
nitudes lower than observed. 

When drift heating corresponding to the velocity of 
only Vd = 2 km is considered, the gas temperature 
becomes much larger than temperatures of the dust. The 
maser optical depth increases by a factor of three and 
the resulting power by two orders of magnitude reaching 
$ = 0.4 cm~^ s~^. This increase occurs because maser 
pumping is the most effective when gas is hotter than the 
dust (see, for example, BP04). 

We can conclude that gas heating by collisions with 
dust (due to thermal motion) is not sufficient to produce 
observed maser luminosity, while additional drift heating 
by moving dust is capable of explaining the masers in AGB 
winds. 



4. Summary 

We have considered the maser effect in a medium con- 
sisting of a mixture of gas (hydrogen and water vapor) 
and dust of various types. The gas and dust temperatures 
and level populations of water molecule are calculated self- 
consistently from the system of population balance equa- 
tions and thermal balance equations for the gas and dust 
in the radiation field of a late-type star. 

When dust of different types is present, the gas inter- 
acting with the grains can be heated by one type of dust 
and is cooled by another. The gas temperature then takes 
an intermediate value. Radiative cooling by water and the 
presence of hot dust strongly influence the water molecule 
energy level populations and therefore should be taken 
into account in calculating the maser effect. We find that 
for a small slab thickness H the inversion appears because 
of the de Jong (1973) mechanism, while for large H, the 
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maser can be pumped by radiation from the dust, whose 
temperature differs from that of the gas. 

The maser strength depends on the combination of 
dust types. When the medium is optically thick to the 
line radiation, the inversion of the 615 — > 623 maser level 
populations appears only for combinations of silicates with 
carbon (or graphite). The main cycle of maser pumping 
is 623 ^^^^ rpj^g upward transition 

in this cycle is dominated by radiation from the hot dust, 
while the heat sink is realized by photon absorption by the 
cold dust. Combinations of water ice with any other dust 
type produces no inversion because of the gap in the ice 
absorption coefficient near 30 /im. Thus, masers operating 
on the differen ce between the gas and ice temperatures 
toeguchilllQSlt BP04) are extinguished by the hot dust 
radiation. 

Strong masers can also be produced if there is a size 
distribution of the dust grains. The maser effect appears 
due to the de Jong mechanism for all discussed dust types, 
if slab half-thickness H < 10^^ cm. We find that the maser 
disappears at iJ < 10^^ cm, because the gas temperature 
becomes too low due to water cooling. For graphite, the 
inversion exists also at 77 > 10^*^ cm, where the maser 
pumping cycle appears to be the same as in the two-dust 
maser model. 

We show that the hot-cold dust model is able to repro- 
duce the strength of water masers observed from a disk 
around the companion of the carbon star in the binary 
system V778 Cyg. However, the masers in the winds of 
AGB stars require an additional source of heating, for ex- 
ample due to friction between drifting dust grains and the 
gas. 
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